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Abstract 

We address the Kondo lattice with one conduction band by introducing a 

representation of the (s = 1/2) local moments consisting in an s = 1 triplet 

of spin fluctuations of a vacuum state. The excitation of these fluctuations 

by the electronic field results in the creation of new fermionic modes. It is 

shown that, in the one-electron case, the zero-temperature spectrum of these 

excitations consists in two s = 1/2 doublet bands and an s = 3/2 quadruplet. 

One of the doublets is seen to be a slight modification of the bare initial 

electronic band for low values of the exchange coupling, while the other bands 

correspond to complex, coherent recombinations of the bare electron and spin 

fluctuations displaying large effective masses. 
PACS: 71.27.+a, 71.28.+d, 75.20.Hr, 75.30.Mb. 
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I. INTRODUCTION 



Heavy-fermion metals display a rich phenomenology produced by the interaction be- 
tween conduction electrons and partially filled /-bands from (typically) Ce or U.Si Its high- 
temperature properties are those of a collection of weakly interacting local moments and con- 
duction electrons with quite ordinary masses. At low temperatures, however, they exhibit 
Fermi liquid behaviour with values of the linear specific heat and temperature-independent 
Pauli spin susceptibility corresponding to huge quasiparticle masses. This Fermi liquid be- 
comes, in many systems, antiferromagnetic, superconducting, or both. 

The high-temperature behaviour described above, suggests that these systems should 
be represented by the Kondo lattice model (KLM), in which the conduction electrons are 
coupled to the local moments through an exchange interaction. This model is the limit of 
the more general Anderson lattice model (ALM) when the fluctuations in the number of 
/-electrons per Ce or U site can be neglected.i 

A number of techniques have been proposed for dealing with these models (especially 
the ALM).B! Of particular importance have been the mean-fieldi and variational^ approx- 
imations. In both cases, the Hamiltonian of the ALM is essentially transformed into an 
effective Hamiltonian in which the /-electrons are uncorrelated, shifted close to the Fermi 
level and hybridize weakly with the conduction electrons. This high density of states at the 
Fermi level would explain the high values of the linear coefficient of the specific heat and 
the magnetic susceptibility. The formulation of these theories relies heavily on the fact that 
the occupation of the /-orbitals is not exactly integer. 

In this article, we present a new approach to the description of the heavy-fermion state. 
By selecting the KLM, we are effectively asserting that the phenomenology of these systems 
is associated with the interaction of conduction electrons and local moments and that the 
/-charge fluctuations do not play an essential role. The general philosophy of our approach 
bears strong parallelisms with the physical discussion of the Kondo model, namely, the inter- 
action of a single magnetic impurity with conduction electrons through an antiferromagnetic 
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exchange coupling.H In this model, the excess magnetic susceptibility is seen to depart, at a 
characteristic (T^) temperature, from the Curie- Weiss law, smoothly approaching a constant 
limit as T — > 0. This is because at T = 0, the local spin forms a collective singlet state 
with the spins of the conduction band with a binding energy of essentially KbT$ Thus, 
for temperatures below Tj,, the resonant scattering of the conduction electrons acts as an 
effective temperature for the local moment giving rise to the mentioned T-dependence of 
the magnetic susceptibility. 

Analogously, in the treatment of the KLM presented in this paper, heavy fermions appear 
as collective, many-particle states made up of electrons and spin fluctuations. Nevertheless, 
while the Kondo resonance is an s = singlet, our collective excitations are fermions. An- 
other difference is that, while the Kondo resonance can form only for antiferromagnetic cou- 
plings, the situation in the KLM is the same for both ferro and antiferromagnetic couplings. 
This is important since heavy-fermion metals like CeAl3 are known to have ferromagnetic 
couplings.0 

The key idea of this article is the introduction of a new, very physical representation 
of the algebra of the local moments. Let us suppose for a moment that the conduction 
electrons are absent. If we consider a lattice of N sites with a local s = 1/2 spin at each 
site, the Hilbert space of the system is spanned by a basis of 2 N states corresponding to all 
the possible spin arrangements constructed from the up and down states of the local spin 
at each lattice site. Since each of these 2 N states has the same energy when conduction 
electrons are absent, the global state of the system can be regarded as a linear combination 
of the 2 N basic vectors, each of them carrying the same probability and an arbitrary, random 
phase. This is what we call the random vacuum state (RVS). We show that by repeatedly 
applying the spin operators y^S^, y2Sj^ and 2S^ to this vacuum state we obtain, in the 
thermodynamic limit, a set of orthonormal vectors that span a new representation of the 
algebra of the local moments. The states of this new Hilbert space can be interpreted as 
spin fluctuations of the RVS. This new representation, initially constructed using physical 
considerations, can be presented at the end in an axiomatic fashion as a perfectly consistent, 

3 



rigorous mathematical structure. 

Thus, with this representation, when conduction electrons are incorporated into the 
system, their interaction with the local moments can be viewed as an electronic excitation 
of spin fluctuations of a vacuum state. In this process, the conduction electrons get dressed 
with spin fluctuations giving rise to new fermionic excitations. For low relative values of 
the exchange coupling (the actual physical situation), we find two types of excitations. One 
type corresponds to slight modifications of the initial bare conduction electrons. The other 
modes are complex, many-particle combinations of electrons and spin fluctuations displaying 
large effective masses. 

We provide a systematic technique to calculate the energy and wave function of these 
excitations. The method essentially consists in determining variationally the best approxi- 
mation to the actual wave function within the space in which the bare electron dresses with 
n spin fluctuations at most. The larger we take n, the better the approximation. 

In this basic, qualitative article, we consider the lowest (n — 1) approximation for the 
system with a single conduction electron only. Although the lowest approximation proves 
to be inadequate for the precise determination of the heavy-fermion bands, it does reveal 
the qualitative picture presented in this paper. The precise calculation of the heavy-fermion 
modes for low relative values of the exchange coupling will require to consider higher orders 
(n > 1) of approximation and promises to be a hard technical task. 

In Section II, we present the detailed construction of the algebra of spin fluctuations. 
The determination of the translational and spin rotational symmetries of the system within 
this new representation is carried out in Section III. The use of these symmetries will be 
essential for the construction and characterization of the fermionic excitations. 

In Section IV we address the simple case of a single conduction electron interacting with 
the lattice of local moments and study how this bare electron gets dressed with spin fluctu- 
ations. The inspection of the atomic limit (band width — > 0) leads to the conclusion that 
the spectrum of the system contains two s = 1/2 doublet bands and an s = 3/2 quadruplet. 
These bands are then calculated in the lowest approximation to determine their qualitative 
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features away from the atomic limit. As we already mentioned, this approximation indicates 
that for small values of the coupling, one of the doublets is a minor modification of the bare 
conduction band, while in the other doublet and in the quadruplet, the conduction electron 
strongly combines with spin fluctuations forming bands with large effective masses. We also 
find in this analysis, that the two doublets hybridize only above a certain critical value of 
the coupling. 

The natural extension of this picture to the case of a system with a finite electronic density 
would require to treat the bare conduction electrons as a Fermi sea with electrons and holes 
as elementary excitations with positive energies. These electrons and holes will be the ones 
giving rise to collective modes when they combine with spin fluctuations. The existence 
of collective states with negative energies would lead to the condensation of some of these 
modes in the ground state, forming the heavy-fermion gas that is observed experimentally. 
The detailed study of the Kondo lattice with a finite density of conduction electrons within 
the framework of the approach introduced in this article is intended to be the focus of our 
future endeavours. 



II. SPIN FLUCTUATIONS OF A RANDOM VACUUM STATE 

We shall study a KLM consisting in a conduction band coupled to a lattice of s = 1/2 
local momemts by an exchange interaction. The Hamiltonian is given by 

H = € k cl a C ka + JS e iSfi , (1) 

where S ei and S/j are the spin of the conduction electrons and the local spin at site % 
respectively. As in many equations throughout this paper, sum over repeated indexes is 
implicitly understood. S ei is given in terms of the electronic operators in the site (Wannier) 
representation by 

^>ei = \c\ a a a pCip. (2) 
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If Rj is the position of the z-site, the c| a -operators are given, for a lattice with N sites, 

by 

4 = iV- 1/2 Ee- iRik cL- (3) 

k 

The fundamental operators c^ a , S/j are characterized by the following algebraic proper- 
ties: 

{eka,Ck'/?} = Skk'SaP , {Cka, Ck'/j} = , (4) 
"S/i> S/j = i8ij£al3'yS'j i , (5) 

S}, = S/, , (6) 
S% = 3/4, (7) 
[cka.S/i] =0. (8) 

As we said before, the fundamental idea of this paper is a novel treatment of the local 
moments: We shall first consider these degrees of freedom only, adding the conduction 
electrons at the end. 

We can satisfy (|5|) and (0) by writing 

S/i = 2^a a a/3fil3 , (9) 

where f} a and /j Q are creation and annihilation operators of /-electrons satisfying cannonical 
anticommutation relations. From Eq. @ we get 

2 3 t 3ft 

= ifiafia + l^fi^filf^fii > (10) 

which implies that Sjj = if the z-site is empty or doubly occupied and S 2 ^ = 3/4 if it 
is singly occupied. Therefore, as it is physically obvious, the condition (|7|) is equivalent to 
demanding that each site be occupied by one and only one /-electron. 

Thus, the Hilbert space in which the S/j-operators act, should be the set of all vectors 
of the form 
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^ Cai,...,ajv AaiJ 2a 2 ■ ■ ■ ZiVajvIO)- (H) 

Since the spin operators S/j preserve the number of /-electrons at each site, they are 
well-defined in this space. It will be convenient for future developments, however, to work 
with the following equivalent set of operators: 

Sl , i = V2(s x fi + iS%)^V2fy il , 

8 0>i = 2S' fi = fy it -flf ili (12) 
8- 1 , i = V2(S%-iS u ft )=y/2flf iV 

The initial S ^-operators can always be recovered using the inverse transformation: 

SXfi = 271 (SM + S ~ m) ' SVfl = 2^2 {Sl ' 1 ~ ' ^ = \ S °' 1 • (13) 



It is straightforward to see that the operators (|i~2) acting on the space satisfy the 
following multiplication table: 



S0,i s l,i 


— Si^ 


s 0,i s -l,i 


s l,i s 0,i 


— — 


s -l,i s 0,i 




= 1 + S ,i 


S 2 




= 1 — s 0,i 


s 2 



(14) 



as well as the properties 



s i,t — s -i,« 5 s o,i — s o,i , (15) 

=0 Vi^ j ■ (16) 

Eqs. (|I~4|)-(|T6|) can also be regarded as the defining conditions of the algebra since, as it 
can be readily checked, they imply Eqs. 

Let us now think in more physical terms: Since the local moments do not interact 
among themselves, all the spin configurations will have the same energy when the conduction 
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electrons are absent. Thus, the physical state of the system can be regarded, in this situation, 
as a random linear combination of the 2 N possible spin configurations, namely, 

\RVS) = ^= £ e«-x-^/J ai ...A aw |0> I (17) 

where the phases R ai ...a N are taken at random. We will call this state the random vacuum 
state (RVS). 

We can now apply to this state our fundamental sjj-operators. As it is obvious from Eq. 



(J12f), their effect is to excite spin fluctuations in the vacuum state. sx,« annihilates all the 
components with spin up at site i and flips the i-spin of the others; the result is a normalized 
state with the z-spin up. s-i^ does exactly the contrary and So,i introduces a relative phase 
shift of 7r between the i-spin up and down components. 

By repeatedly applying the operators S[j to the RVS, we generate new states like, for 
example, 

|V>> = s 1)2 s 0)5 s 1)15 \RVS). (18) 

We shall loosely say that has spin fluctuations at sites 2, 5, and 15. When applying 
Si t i to a state that already has in its expression an operator in the same lattice site, we 
can use the multiplication table (|14|) to produce a state with one s^-operator at most. For 
example, 

S-l,l5\i>) = S lj2 SQ t5 \RVS) - Si t2 So,5So,15\RVS). (19) 

We can construct 4^ different combinations of spin fluctuations by leaving each site 
either empty or with one of the three possible fluctuations. The linear combinations of these 
states form a vector space in which the s/j-operators are perfectly defined through Eqs. (|I4]) 
and (0). 

Using the randomness of the vacuum state, it is straightforward to prove that the 4^ 
basic states constructed above are orthonormal in the thermodynamic limit. For example, 
the scalar product of so^\RVS) and s^i\RVS) is given by 
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1 1 



E e 




(20) 



V2 2 N -i 



a 2 ,...,ojv=T,-L 



which is essentially the average of 2 random phases and, therefore, vanishes as iV — > oo. 
The result for the general case can be obtained in an analogous way. 

It should be emphasized that the exact orthonormality holds only in the thermodynamic 
limit (N — ► oo). Actually, for N finite, it would be mathematically impossible to have a set 
of 4^ exactly orthonormal states in a vector space of dimension 2 N . 

At this point, the physical ideas introduced in this section and in particular the concept 
of RVS can be regarded as guiding physical tools in the construction of our new algebraic 
structure. Now that this structure has been established, we can completely ignore the 
concept of RVS and present, as we do next, the obtained representation of the algebra of 
the local moments in an axiomatic, mathematically rigorous way: 

1. Vector space: It is the set of all states of the form 



where C hh ...i N G C, u 2 ,i = 1, and u hji = s hji for k = -1, 0, 1. 
2. Operators sj j : Their action on the vector space is given by the multiplication table 



3. Scalar product: It is defined by declaring that the basis of the space formed by the 4 
vectors i% i ittz 2i 2 • • • Ui n ,n\&) is orthonormal. 

These three points define a Hilbert space with operators satisfying Eqs. 0l4[)-(|l6|) or, 
equivalently, with operators S/j [given by Eq. ([13])] satisfying the defining conditions of the 



The inclusion of the conduction electrons is straightforward: The basis for the complete 
Hilbert space is obtained by applying electronic creation operators to the basis of spin 
fluctuations presented above, and the action of the operators cj^, Ck« on the states is given 



E c, 



'l 1 l 2 ...l N U ht iU ht2 ■ • .Ul mN \$), 



(21) 



ii=-l,0,l,2 



(UP and Eq. (fT|). 



algebra 
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by the cannonical anticommutation relations (f|), Eq. (H) (where we can substitute S/j by 
s^i), and the condition Ck a |$) = Vk, a. 

Finally, the Hamiltonian is given, in terms of the basic operators of this new representa- 
tion, by: 



J 



(22) 



III. TRANSLATION AL AND SPIN ROTATIONAL SYMMETRIES 

In this section, we determine the translational and spin rotational symmetries of the 
system within this new representation. These symmetries will be used in the following 
sections to characterize the spectrum of excitations of the system. A unitary operator is 
completely determined by specifying how it transforms the basic operators c^ a , s^i and its 
action on the vacuum state: this is what we do next to introduce the symmetries. 



A. Lattice translations 

It is physically evident that, if a 7 (7 = 1,2,3) are three primitive vectors of the real 
lattice, their associated translation operators T aj should satisfy: 

^a 7 CR ! Q T^ = c R . +a ^ Q , (23) 
T^si^T^ = s ZjRi+a7 , (24) 
T !h \$) = m. (25) 

These equations define unitary operators that leave the Hamiltonian invariant and thus 
represent the basic translational symmetries of the system. 

To exploit this symmetry it is convenient to Fourier transform all operators to the mo- 
mentum representation: 

cL = A^ 2 £e^ k cL (26) 

i 

Sim = N- 1 ' 2 £ e^ k s M , (27) 
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where k is always in the first Brillouin zone. Lattice translations just multiply these operators 
by a phase factor: 

T a A«Tl = , (28) 

T ai s lM Tl = e-^s^ . (29) 

The s/ k operators satisfy the following equations, which will be of much use in the next 
sections: 

[so,k, s±i,k'] = ±2iV- 1 / 2 s ±lik+k , , (30) 
Kk, s-i,k'] = 2iV- 1 / 2 s 0ik+k , , (31) 
kk,s/,k'] =0, (32) 
4k = s -i,-k • (33) 



Finally, the Hamiltonian is given, in the momentum representation, by 
J 

AN 1 / 2 



H = £kC ka c ka + ^ /2 [(c kT C k / T - C ki C k /|)s , k '- k + V^cj cT CVlS-l,k'-k + V^C^Ck^S^-k 



(34) 



B. Spin rotations 



The spin rotational symmetry is quite subtle. The total-spin operators of the system are 
given by 



S — S ei + Sfi , 



(35) 



and satisfy 



s f = s, 

[#,£] =0. 



(36) 
(37) 
(38) 
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These are the properties that characterize the generators of a unitary representation of 
SU(2) that leaves the Hamiltonian invariant. There is, however, an important problem: the 
operators S do not annihilate the vacuum state and, therefore, the unitary transformations 
generated by them would not leave |<&) invariant. Since (in the absence of magnetic order) 
we do not physically expect the spin rotational symmetry to be broken, X cannot be taken 
as the generators of this symmetry. We should thus look for another set of operators S 
satisfying Eqs. (|36|)-(|3~8l) plus the condition 

S|$) = 0. (39) 

As we already saw, every vector in the Hilbert space constructed in Section II can be 
written as V|$), where V is a linear combination of products of operators c\ a and sy. The 
effect of a transformation exp(zaS) on a state would be 

e iaS V|$) = e iaS Ve- iaJ: [e ioS |$)] , (40) 

where 

e iQS |$) |$). (41) 

The obvious thing to do to satisfy Eqs. (j36|)-(|38|) plus the condition (j39"D, is to demand 
that exp(iaS) transform V in the same way as exp(iaS) but leave |$) invariant, namely, 

e iaS V\$) = e iaS Ke" iaS |$), (42) 

which implies that the action of S on the Hilbert space is given by 

&V\Q) = [S,F] |$). (43) 



It is quite straightforward to prove that the linear operators S defined by Eq. ( [43]) satisfy 
Eqs. (|36|)-(|3"9|) and will, therefore, be taken as the generators of the spin rotational symmetry 
of the system. This implies, in particular, that the conduction electrons 0^1$) (a =t, I) 
and the spin fluctuations sj,k|$) (I = — 1, 0, 1) form, respectively, an s = 1/2 doublet and an 
s — 1 triplet under spin rotations, as is physically expected. 
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IV. KONDO LATTICE WITH ONE CONDUCTION ELECTRON 



We shall consider, in this section, the case of a single conduction electron interacting 
with the lattice of local moments. This system is interesting because it is simpler than the 
general case and reveals very clearly the appearance of heavy fermions as collective states 
of the bare conduction electron and spin fluctuations. 

We shall first examine the atomic limit (band widths 0) to determine that the spectrum 
of the system consists in two s — 1/2 doublet bands and an s = 3/2 quadruplet. Then, 
we will study the structure of these bands away from the atomic limit using the lowest 
approximation in which the conduction electron can only dress with one spin fluctuation at 
most. Finally, we will see that this study reveals the formation of collective, heavy- fermion 
states for low relative values of the exchange coupling. 

The KLM is particularly simple in the atomic limit; it is just the sum of local Hamilto- 
nians associated to each lattice site: 

H = Y J H l (Hi = JS ei S fi ). (44) 

i 

Thus, in this limit, the one-electron problem reduces to diagonalizing Hi in the space of 
the local degrees of freedom with a single electron present. In the representation introduced 
in Section II, this space is generated by the following eight vectors: 

a=U, (45) 
4 SM |$) a=U, 1 = -1,0,1; (46) 

(a state with more than one spin fluctuation accompanying the conduction electron at site % 
can be developed by the multiplication table (|14"D into a linear combination of these states). 

Since the spin operators S are well-defined in this local space and commute with Hi, they 
can be used to diagonalize the local Hamiltonian. The diagonalization of S 2 is straightfor- 
ward: the vectors ([45]) form an s = 1/2 doublet while the six vectors (|46| ) are the product 



of an s = 1/2 by an s = 1 representations and can be decomposed into an s = 1/2 doublet 
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and an s = 3/2 quadruplet. The energy of the quadruplet is seen to be J/4 while the 
diagonalization of if; in the subspace of s = 1/2 leads to two doublets with energies J/4 
and —3 J/4. The states in each multiplet can be determined, for example, by repeatedly 
applying the lowering operator (S~ = S x — iS y ) to the states with the highest eigenvalue of 
S z (s z = s). The explicit expressions of these vectors are given below. 

1. Quadruplet (E = J/4): 

\j,l,l,i) = 4, Sl M, (47) 
= ^-\ J -\, | = ±(cl lSl ,.-V2cy , m , (48) 



J 3 -1 
4' 2'T" 



,i) = -\ S ~\j>l> = 7f (c * f T S -M + ^4^0,)!$), (49) 

iff'T'^^lff'T'^^-^- (50) 



2. E = J/4 doublet: 



'l' ^ ^ = IT (3 ^ T + ^ + ^Vm)!*), (51) 
if ^ ^.0 = ^f ^ ^> = - <ti*o,i + V2cJ T *-i,i)l*>- (52) 



3. E = -3J/4 doublet: 



-3J 1 1 



11 1 

l^'^^pO = ^l-^'^'*) = \( C K + c U s o,i - v^ct T s_ M )|^>. (54) 

Thus, the Fourier transform of these vectors, namely, 

|.E7, s, s 2 , k) = N~ l l 2 e iR *\E, s, s z , i) (55) 

i 

(with k in the first Brillouin zone), are exact eigenstates of the operators if, S 2 , S z , and T &7 , 
in the atomic limit, with eigenvalues E, s(s + 1), s z , and exp(— ika^) respectively. These 
operators form a complete set in the sense that each vector is characterized by a definite 

14 



set of eigenvalues. Since S 2 , S z , and T ai always commute with H, they can be used to 
characterize the spectrum of the system away from the atomic limit as well. This simplifies 
enormously the calculation of the spectrum: we just have to diagonalize H in the subspace 
of vectors with definite eigenvalues of these operators. 

Thus, from the analysis of the atomic limit and the symmetries of the system, we conclude 
that the spectrum contains two s = 1/2 doublet bands and an s = 3/2 quadruplet. We shall 
now propose how to calculate these bands away from the atomic limit. Since all the members 
in a spin multiplet have the same energy, we only need to calculate the energy of one of 
them; we shall select the one with s z = s. 

The most general state of the one-electron system is a linear combination of terms con- 
taining a single cj^-operator multiplied by several spin fluctuations corresponding to various 
lattice sites. The approximation that we propose to calculate the spectrum of excitations 
consists in determining variationally the best approximation to the actual eigenstates of H 
within the subspace in which only a maximum of n spin fluctuations can accompany the 
conduction electron. 

One should start by writing the most general vector in this subspace with definite eigen- 
values of T &7 , S 2 , and S z . We shall denote this state (with s z = s) by |£, s, s, k) = V£ )S) k|3>), 
where £ represents the parameters of this general vector. The diagonalization of T &7 is 
straightforward if we consider all the degrees of freedom in momentum space, and the spin 
eigenstate conditions are equivalent to 



where | _L) is a vector orthogonal to |£,s,s,k). Obviously, the best approximation to the 
actual eigenstate will be the normalized vector for which (_L | _L) is minimized. Thus, the 
optimal ^-parameters should be determined by demanding that 




(56) 



S z \£,s,s,k) = [i; z ,Vs, s , k ]m = s\Z,s,s,k). 



(57) 



The action of H on |£, s, s, k) can be written as 




(58) 
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(±\±-) = (Z,s,s,k\H%s,s,k)-El (59) 
is minimized, where _E k (the approximate energy) is given by 

£ k = (£,s,s,k|iJ|£,s,s,k), (60) 
and |£, s, s, k) satisfies the normalization condition 

<£, S)S ,k|£, S , S ,k) = l. (61) 

In this article, we will consider only the lowest (n = 1) approximation. From Eqs. (f4T|)- 
(|55|) it is clear that this approximation gives the exact result in the atomic limit and it should, 
therefore, yield precise results around this limit. It also describes exactly the bare electrons 
when J = and it is thus expected to express accurately how these excitations are modified 
when a small interaction is introduced. Away from these situations, this approximation is 
expected to describe the spectrum only qualitatively. 

The determination of the states |£, 1/2, 1/2, k) and |£, 3/2, 3/2, k) in the lowest approx- 
imation is straightforward, the result being 

\[a,0\,^, ~ k) = [« k cj cT + iV- 1 / 2 /? k (q)( c ; TS o,k- q + v^c^k-q)] |$>, (62) 

|[7]4| k ) = Ar_1/2 ^(q)4^l $ )' ( 63 ) 

where sum over q is implicitly understood. It is also assumed that k — q is in the first 
Brillouin zone, which occasionally requires the addition of a reciprocal lattice vector, a, 
(3, and 7 are the variational parameters whose optimal values should be determined by 
the conditions (j59"|)-(|6~l|). The mathematical development of these conditions is a lengthy 
technical task which is outlined in the Appendix. 

Since we will be considering the system in all the situations between the atomic and 
J = limits, it is important to write the Hamiltonian in terms of parameters that scale 
in a convenient way. Since the number of conduction electrons is a constant, we can shift 
the bare dispersion relation so that = Ar] k , where A is the band width and r/ k is a 
normalized band of unit width centered at e = (—l/2<r]] i < 1/2). 
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Instead of working with J, A, and r^k, we will use the more convenient parameters A, A, 
and ?7k given by 

x -jdwr (64) 

A = A + |J|. (65) 

|A| is a dimensionless parameter describing the strength of the interaction and it varies 
from zero (no interaction) to one (atomic limit). Positive and negative values of A correspond, 
respectively, to antiferromagnetic and ferromagnetic couplings. A is just a measure of the 
global energy scale of the system which factorizes in H, and can thus be set equal to one 
without loss of generality. The Hamiltonian, therefore, depends essentially on the parameters 
A and r/ k and is given by 

H = e k cL C ka + ASeiS/j, (66) 

where e k = (1 - |A|)r/ k , -1/2 < r? k < 1/2, and -1 < A < 1. 

As we said before, the development of the conditions (|59|)-(|6lD that determine the optimal 
parameters of the approximate eigenstates (|62"D and (|63f) is carried out in the Appendix. We 
next summarize the results obtained for the doublet and quadruplet bands: 



1. Doublet bands: The optimal parameters a and (3 in Eq. (|62|) and the approximate 
energy E are given by 

a k = Re [4a k A k ] , (67) 

Ac(q) = X \ + ^ ; , (68) 

E k = e k - Re[w k ] , (69) 



where 



of = L . i (70) 

wA Im[w*A*[8Aw*w-(4:A + 3B/A)w*(w*-w)-6\w)) ' 

A>^r n ° Md \ -2. (71) 

B = A 2 C' 2 D{v)dv (72) 

k l 1/2 [(l_|A|)7 ? -£ k + U ; k ] 2 ' 1 J 
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and Wk is the solution of the following pair of equations 



Im 



w*(8w*w + 3A(1 + 2/A)w* - 3X 2 /2) 



(73) 



Im 



w* 2 A* 



(AA + 3B/A){w* - wf + 3A(A + 2){w* - 3w) + 6A(1 - A/A*)w 



-- 0. 
(74) 



D{rf) in Eqs. ( [fTD and ( \T2\) represents the density of states of the normalized band 
77k- The subscripts k in Eqs. (|70|) , (|73|) , and (|74]) have been omitted to simplify the 
notation. 

2. Quadruplet band: The optimal 7 in Eq. (^) and the approximate energy E are given 
by 



7(q) 



a a 

£q Z £q Z 



E = Re{z] 



where 



a 



\C*lm[z] 



2CIm [C*(l + Dz/C)\ 



-1/2 (1 - |A|)»7-z 



L> = A . 

-1/2 [(l-|A|)i7-s. 
and z is the simultaneous solution of 



2 ! 



Im 



, Im [C{zC* - 3A)] = 0. 



(75) 



(76) 
(77) 
(78) 



(79) 



It is clear from these equations that the energy E^ of the doublet bands depends on A, 
£k, and D{rj). In the case of the quadruplet, however, E^ depends exclusively on A and 
D(rj), which implies that this band is flat in this approximation. 

In order to carry out explicit calculations, we must specify D(rj). A realistic density of 



states should behave like J 1/2 ± 77 at the extreme points 77 = =Fl/2 and must be normalized. 
We shall take the simplest function satisfying these conditions, namely, 



D{ V ) = ^sj\- V \ (80) 



which leads to simple expressions for A, B, C, and D [see Eqs. ( |A34 ) and ([A3 5| ) ] . 



At this point, we can already determine the lowest approximation to the spectrum of the 
system by numerically solving Eqs. (|73|), flT4"|) and ([T9|). 

The solutions for the energies satisfy E(— A, — e) = —E(\,e) for the doublet bands, and 
E(—X) = —E(X) for the quadruplet. Since we can relate, with these equations, the spectra 
for positive and negative values of A, we will restrict our study to the antiferromagnetic case 
(A > 0). 

The numerical results obtained are displayed in Fig. 1. In order to make the physical 
interpretation of the results easier, we have parameterized the values of the bare band e in 
a band-like manner, namely, with the function e(9) = —1/2(1 — A)cos# [9 G [0, 7r]), which 
is represented by the dashed curve. The solid curves correspond to the solutions obtained 
for the two doublets E(\,e(9)) and the quadruplet E(X). 

Fig. 1(a) represents the atomic limit, with completely localized bands at E = 1/4 
(quadruplet and upper doublet) and E = —3/4 (lower doublet). In Fig. 1(b), the degeneracy 
of the upper bands in the atomic limit is broken; as we already mentioned, the quadruplet 
band remains flat for all values of A in this approximation. The transition from Fig. 1 (c) to 
Fig. 1(d) is very revealing: There is a clear tendency to form, with the states of the superior 
and inferior parts of the upper and lower doublets respectively, a band that closely follows 
the bare conduction band, while the rest of the states tend to form a practically flat band. 
We shall call these bands the continuous (because it is essentially a continuation of the bare 
band) and the heavy-fermion band respectively. Fig. 1(d) can actually be interpreted as a 
hybridization of these two bands with the characteristic pseudogap structure. 

To get a deeper insight into the nature of these bands, we will study the parameter a. It 
is clear from Eq. fl62|), that a£ represents the probability for the mode | [a, /?], 1/2, 1/2, k) to 
be in the bare state cL|$). If the subscripts U and L denote the upper and lower doublets 
respectively, we get, for A = 0.26, 
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a\ (9 = 0) = 0.77 , al{6 = vr) = 0.93, (81) 
al(6 = 0) = 0.22 , a 2 L (9 = tt) = 0.018. (82) 

The points in rows and (|82| ) correspond to the continuous and heavy-fermion bands 
respectively. The interpretation is clear: While the continuous band consists essentially in a 
small modification of the bare electronic states, in the heavy-fermion modes, the conduction 
electron strongly combines with spin fluctuations forming a complex, collective state as in 
the Kondo problem. It should be noted that the s = 3/2 quadruplet states are also collective 
modes in the most strict sense since, due to spin conservation, they have no probability to 
be in the bare states. 

Calculations for very small values of A show that the continuous band actually coincides 
with the bare band as A — > 0. Since, as we already argued, these calculations are very reliable, 
we can assert that there is no hybridization of the continuous and heavy-fermion bands for 
very small couplings. Figs. 1(f), 1(e), and 1(d), describe the process of hybridization of these 
bands around A = 0.24. We find that beyond the poins A, B, C, and D in Figs. 1(e) and 
1(f), the system of Eqs. (|73|) and (|74]) no longer has a solution. We have interpreted this fact 
as a mathematical expression of the inadequacy of the lowest order approximation around 
these points: one should physically expect that the precise description of the collective states 
will require considering more than just a single spin fluctuation combining with the bare 
conduction electron. 

V. SUMMARY AND CONCLUSIONS 

With the introduction of a new representation of the algebra of the local moments, we 
have established a new framework for the study of the Kondo lattice which enables a very 
clear mathematical description of heavy fermions as collective states of conduction electrons 
and spin fluctuations of the lattice of local moments. 
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APPENDIX: LOWEST ORDER EQUATIONS FOR THE DOUBLET BANDS 



We shall briefly develop here the conditions (|59|)-(|6lD for determining the optimal param- 
eters a and (3 that characterize the doublet bands of the one-electron Kondo lattice in the 
lowest approximation. The development of the corresponding equations for the quadruplet 
band is completely analogous and will be omitted. 

Let us summarize the problem: We have to determine the parameters a and (3 of the 
vector | [a,/?], 1/2, 1/2, k) in Eq. fl62|) that minimize the function 



* k ([a, fl) = ([a, 0\, i, i, k\H 2 \ [a, /3], i, 1, k) - E* , (Al) 

where 

E k =([a,/3],^,k|F|[«,/3],i,i,k), (A2) 

while satisfying the normalization condition 

([a,/3],i i k|[a,/?],i 1 k> = l. (A3) 

In order to handle this normalization condition, it is very convenient to introduce a 

Lagrange multiplier r. Thus, the objective will be to find the parameters a and (3 that 



minimize 



n k ([a,(3],T)=^([a,(3})+T([a,(3},± i k|[a, /?], ~ ~ k), (A4) 



for an arbitrary r, and then select the value of r for which the normalization condition ( |A3| ) 
is satisfied. 
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To simplify the notation, we will often omit the subindexes k in the forthcoming equa- 
tions. If mi, m 2 , ni, and n 2 are defined by 



(A5) 



mi = f //3(q)dq, m 2 = J /^*(q)/3(q) dq, 
m = i/e q /3(q)dq, n 2 = J / £ q /3*(q)^(q) dq , 

where the integrals are extended over the first Brillouin zone and T is the volume of this 
zone, a lengthy calculation yields 



11 11 

([«>£]> 2' 2' k l[ a '^' 2' 2'^ = a * a + 3m2 ' 



(A6) 



, 3A , . 3A . 

E = a ae H (anij + a m 1 m 1 + 6n 2 

4 2 



(A7) 



-5 f iAr? 
([a, /?], -, -, k|tf 2 | [a, /?], -, -, k)= - J e\{? (q)/%) dq + -^(a/2 - m x ) + 

3An! 



9 3 

fa*/2 - m?) H A 2 m 2 + -A 2 mtmi + 

2 v / i; 16 2 4 1 

3 

~A(e - A/2)(am£ + a*mi) + (e 2 + 3A 2 /16)a*a . 



(A8) 

Substituting these equations in (|Aj), we can determine fl^ ([a, (3], r) and calculate its 
variations with respect to /3(q) and a. The equations that express the vanishing of these 
variations are, respectively, 

3A 2 



= /3(q 
Xa , 



<-2Ee q + — + r 
Ani 



A . A 2 m x 

+ -(a/2-m 1 )£ q +^p + 



{ x e-X/2)-^-XE(a/2- mi ) Vq, 

3A 

= a (e 2 + 3A 2 /16 + r - 2Ee) + — [m + mi(e - A/2 - 2£)] . 



(A9) 
(A10) 



We see from Eq. (|A9|) , that the function /3(q) has the structure 



Xa Xb 
+ 



£q % £ q Z 



(All) 



where the complex parameters a, b, and z satisfy 
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a + b=-- + ^-, (A12) 

. zb - z * a= r h.- ?j±( z + z * + a/2) + j( z + z* + A/2 - e) , (A13) 

3A 

z + z* = 2E = 2aa*e + — (ami + a* mi) - 3Am^mi + 6n 2 , (A14) 
lb 



A complete set of equations is formed by Eqs. ( |A12|) - (|A15| ) plus Eq. ( |A10| ), which can 
be written as 

a{e -z)(e- z*) + -j [m - [z + z* + A/2 - e)} = 0, (A16) 
plus the normalization condition 

a*a + 3m 2 = 1. (A17) 

Manipulating these equations and reparameterizing them with w = e — z, we get the 
following equivalent set of equations: 

= a(2w-\)+b(2w*-\)-mie + m, (A18) 

3A 3A 
= — [ a (w - A/2) + 6(V - A/2)] + 4w*w(a + b- m^/2) - —(w + w* - X/2)m u (A19) 

= w * + w - 6{em 2 - n 2 ) - 3A [2{a* + 6*)mi + 2m*(a + 6) - mj mi ] , (A20) 

= ^ - 4(a* + 6*) (a + 6) - ^ - rn^rni + 2 [(a* + &*)mi + m\(a + 6)] , (A21) 

a = 2mi -4(a + 6), (A22) 

r = ( £ -«;*)(£-«;)- — . (A23) 

The first four equations of this set should determine the parameters a, b, and w, while 
the last two equations simply give us the values of a and r in terms of these parameters. 
To get closed relations between a, b, and w, we have to substitute in Eqs. ( |A18|) - (|A21|) the 
explicit expressions of mi, m 2 , ni, and n 2 in terms of a, 6, and u>, namely, 

m 1 = aM + bM*, (A24) 
m = a [A + (e - w)M] + b [A + (e - w*)M*] , (A25) 
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A- 

w* — w 

where M and B are given by 



m 2 = a *bB* + VaB + A a * a + 5 * 6 ( M _ M *\ 

w* — w 

n 2 = ab* [AM + (e - w)B] + a*b [AM* + (e - w*)B*] + 
a* a + b*b 



[(e - iw)M - (e - w*)M*] 



(A26) 



(A27) 



r7 e 



e + w 



dq 



r J (e q - e + to) 5 



(A28) 



The result of this substitution is, after some manipulations, 



= Im 
= Im 



w*(8w*w + 3A(1 + 2/A)w* - 3A 2 /2) 



(A29) 



a*a 



W * 2 A*[(AA + 3S/A)(tw* - wf + 3A(A + 2)(iu* - 3w) + 6A(1 - A/A*H] > ( A3 °) 
Im[A1 (A31) 

(A32) 



Im [wM*(8Au;*t/; - (4A + 3B/A)w*{w* - w) - 6Xw)} ' 
wA 



w*A* 
where A = M - 2. 

It should be noted that Eqs. (|A29|) - ([A32| ) determine everything except the phase of a. 
This is understandable since a change of this phase leads to the multiplication of the state 
by a global phase factor. Thus, we can fix the phase of a by demanding, for example, that 
b = a*. This selection yields a real wave function and gives rise to the final equations listed 
in Section IV. 

A and B are given, in terms of the density of states D(j]) of the normalized band rj^, by 
"i/2 D(r])dr] 2 ri/2 D(rj) drj 



A = A 



1/2 (1 — |A|)?7 — e + w 



— 

-1/2 (1 



/2 [(1 - |A|)7j-e + w]' 



(A33) 



For the model density selected in Section IV (D{rj) = 1/4 — rf), these expressions 
can be integrated analytically. We get, for Im[w] < 0, 

8A 



.4 
B 



(1-|A|) 2 

-8A 2 
fl-IAI) 2 



w-e + iyj{l- \X\) 2 /A-(e-wf 



i[e — w) 



^(i-\x\y/A-(e-wy 



(A34) 
(A35) 
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FIGURES 

FIG. 1. Zero-temperature spectrum of the antiferromagnetic Kondo lattice with one conduction 
electron, in the lowest approximation, for decreasing values of A. The dashed curve parameterizes 
the values of the bare conduction band and the solid curves correspond to the quadruplet (horizontal 
line) and the two doublets. For low values of A, the doublets give rise to a band that closely follows 
the bare conduction states and a practically flat, heavy-fermion band. 
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